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Abstract 

In the atmospheric neutrino experiments the primary problems are the huge uncertainties of flux, very rapid 
fall of flux with increase of energy, the energy dependent wide resolutions of energy and zenith angle be- 
tween true neutrinos and reconstructed neutrinos. These all in together make the choice of binning of the 
data for chi-square analysis complicated. The large iron calorimeter has the ability to measure the energy 
and the direction of the muon with high resolution. From the bending of the track in the magnetic field it 
can also distinguish its charge. We have analyzed the atmospheric neutrino oscillation generating events 
by Nuance and then considering the muons produced in the charge current interactions as the reconstructed 
neutrinos. This practically takes into account the major problem of wide resolutions. We have binned the 
data in three ways: i) in the grids of \ogE — logL plane, ii) in the grids of \ogE — cos 9zenith plane, and iii) 
in the bins of log(L/£'). We have performed a marginalized study over Amgj, 613 and 823 for neutrinos 
and anti-neutrinos separately for each method and finally compared the results. 
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1 Introduction 

The atmospheric neutrino anomaly was first observed by 1MB in 1986 and then confirmed by Kamiokande in 
1988 mill. Finally, the neutrino oscillation was discovered in 1998 with atmospheric neutrino experiment ||3|. 
The atmospheric neutrinos are produced by the interactions of the cosmic rays with the atmosphere. At the 
neutrino energies above a few GeV, the effect of geo-magnetic field on the cosmic rays is negligible and then 
the atmospheric neutrino flux can be predicted to be up-down symmetric. The flight lengths for up and down 
going neutrinos are very different. The atmospheric neutrino experiments exploit these features to study the 
neutrino oscillation. 
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The atmospheric neutrinos are also equally important in the precision era of neutrino physics. The main thrust is 
now on the precise measurements of oscillation parameters. This helps to identify the right track to understand 
the underlying principle that gives the neutrino masses and their mixing. In the recent years, the studies of 
neutrinos has become a popular tool to probe the physics beyond the standard model. In the standard oscillation 
picture, there are six parameters. The present la, 2o and 3a confidence level ranges from global 3v oscillation 
analysis (2008) i [4] are tabulated in table[I] 

Parameter Am^j/lO"^ eV^ sin2ei2 sin^Oia sm^Q23 |Amfi |/10"^ eV^ 

Best fit 7.67 0.312 0.016 0.466 2.39 

la range 7.48-7.83 0.294-0.331 0.006-0.026 0.408-0.539 2.31 -2.50 

2a range 7.31 -8.01 0.278-0.352 < 0.036 0.366-0.602 2.19-2.66 

3a range 7.14-8.19 0.263-0.375 < 0.046 0.331-0.644 2.06-2.81 



Table 1: Global 3v oscillation analysis (2008) 



This spectacular achievement is very stimulating to uncover the facts which are still missing. To determine the 
mass ordering (sign of Am32)B the values of 613 and 8cp with good precision, the octant of 623 with atmospheric 
neutrinos as well as neutrinos from artificial beams, there are many ongoing and planned experiments: INO 
Q, UNO in, T2K Q, NOvA m, Hyper-Kamiokande S and many others. In the current few years, a large 
fraction of effort in particle physics research has gone to study the physics potential of these detectors |[T0| . The 
current research activity lfTTl[T2l[T3l[T4l[T5l[T6llT7l[T8l[T9ll20l shows the uniqueness in physics potential of the 
large magnetized Iron CALorimeter (ICAL) detector at the India-based Neutrino Observatory (INO). It should 
be noted here that its position at PUSHEP has a special feature. It gives the magic baseline from CERN for 
beam experiments, which provides the oscillation probabilities relatively insensitive to the yet unconstrained 
CP phase compared to all other baselines and permits to make the precise measurements of the masses and their 
mixing avoiding the degeneracy issues [16 1. On the other hand, ICAL can detect and separately using the 
magnetic field for charge current events. The oscillation study with atmospheric neutrinos is the primary goal 
of ICAL at INO. Before going into the detailed techniques of the analysis methods, we will first discuss the 
basic nature of atmospheric neutrino oscillation and the detection characteristics of ICAL detector. 



1.1 The atmospheric neutrino oscillation and the ICAL detector 

The present atmospheric neutrino data from the pioneering Super Kamiokande (SK) experiment are well ex- 
plained by two flavor oscillation ||2T1|221 . However, one expects the reflection of — > oscillation in data for 
standard 3-flavor framework in the data if 613 is nonzero. Neglecting the Aniji term the oscillation probability 
can be expressed as: 

P(v„^v,) = P(v,^v^) 

^The CP-violating pliase Sep is still unconstrained. 

■^Am^2 = — 



2 



sin^ 623 sin^ 2613 sin^ 
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X sin^ I \ 



These oscillation probabilities are derived for vacuum. Since it involves electron neutrino, the oscillation will 
be modulated by the matter effect Il23ll24l . Then, 



l.nAm^ML 



P(V^^V,) = P(V,^V^ 



sin" 023 sin" 2Q'l\ sin" ( ^ ) . (2) 



The symbol 'M' denotes effective parameters in matter. The effective mixing angle is 

sin^ 2013 

(cos 2013 -Acc/Am2)2 + sin2 20i3 

and 



sin2 20f3 = "-^^^^^^ (3) 



Arn^M = V (Am^ COS 2013 —Acc )^ + (Am2sin20i3)^ (4) 
with 

Acc = iVlGpNeE, (5) 



where Gf is the Fermi constant, Ne is the electron density of the medium and E is neutrino energy 11251 . The 
matter potential term Acc has the same absolute value, but opposite sign for neutrinos and anti-neutrinos. The 
Mikheyev-Smimov-Wolfenstein (MSW) resonance occurs when neutrino passes through the matter (see eq. (3]). 
It happens for Normal Hierarchy (NH) with neutrinos and for Inverted Hierarchy (IH) with anti-neutrinos. The 
resonance energy corresponding to a baseline can be seen in |26]. 

The muon neutrino (anti-neutrino) produces fi^ (/j+) in Charge Current (CC) weak interactions. The magnetized 
ICAL can distinguish ^+ and /u^ with the magnetic field. The energy (E) and zenith angle (0zenith) or baseline 
(L) resolutions of the muons are very high at ICAL f5\. The hadron energy can also be measured at ICAL. 
However, its resolution is very poor and strongly depends on thickness of the iron layers. 

The atmospheric neutrinos are expected to be very useful in precision studies for its very wide energy range 
(MeV — few hundred GeV) and wide baseline range (few km —12950 km). It gives both neutrino and anti- 
neutrino, which behave oppositely with matter. This helps to detect the sign( Amjj), the value of 0i3 as well 
as the octant of 023. One can exploit this feature to measure the precision of these parameters and the mass 
ordering at the magnetized ICAL detector at INO. It should be noted here that the non-magnetized detectors, like 
water Cherenkov detector, can also contribute in this study since the cross section, the y{= {E^J — £'iepton)/£^v) 
dependence of the cross section are different for v and V. The water detectors may also be able to distinguish 
statistically and due to different capture rates and lifetimes of the charged muons in water. 
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However, one of the crucial problems in neutrino physics experiments is the wide resolutions of E and L be- 
tween true neutrinos and reconstructed neutrinos, which smears the oscillation effect to some significant extent. 
This arises mainly due to interaction kinematics. The un-observable product particles, un-measurable momen- 
tum of recoiled nucleus are the main sources of this huge uncertainty in reconstructed neutrino momentum. 
These are strongly neutrino energy dependent. 

Due to the above complications, the method of extraction of the results from the data is not straightforward. 
The results depend crucially on the way of the analysis and particularly on the type of binning of the data. This 
fact is well-known from the analysis of atmospheric neutrino data of SK experiment 11211 122|| . In this paper we 
consider the reconstructed energy and the direction of an event only from the muon generating it by the neutrino 
event generator Nuance-v3 |[27l . The addition of hadrons to the muon, which might increase the reconstructed 
neutrino energy resolution, is not considered here for conservative estimation of the sensitivity. It would be 
realistic in case of GEANT-based studies since the number of hits produced by the hadron shower strongly 
depends on the iron thickness. However, INO can also detect the neutral current events. Though it is expected 
that these will not have any directional information, the energy dependency of the averaged oscillation over all 
directions can also contribute to the total separately in the sensitivity studies. Here we have studied the 
atmospheric neutrino oscillation by binning the data in different ways and finally compared the results. These 
are discussed in the next sections. 



2 The X analysis 

Now we will describe a general expression for x^, the method for generation of the theoretical data, the esti- 
mated systematic uncertainties, and finally the ways of binning of the data. The number of events falls very 
rapidly with the increase of energy and the statistics is very poor at high energy. However, the contribution to 
the sensitivities of the oscillation parameters is significant from these high energy events. To incorporate these 
events at high energy, the value is calculated according to Poisson probability distribution. For all types of 
binning, we define a general expression of x^ as 



N 



1=1 



= I 2{A^f-A^;}-2A^;in 



Nl 
Nf 



k=l 
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(7) 
(8) 



The N"j (Nfj) is considered as the number of observed (predicted) events in the ijth grid in the plane of logE 



cosSzenith- Here we consider the data for 1 Mton.year exposure of the detector. The f/ - is the systematic error 
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of Nfj due to the kth uncertainty. The is the pull variable for the Mi systematic error. We consider = 5. 
Here we have considered 30 bins of logfi and 300 bins of cosOzenith for both Nfj and Nfj. However, it should be 
noted here that in calculation of the oscillated flux we consider 200 bins of logE and 300 bins of cosOzenith to 
find the accurate oscillation pattern. We consider the E range 0.8 — 50 GeV and cos Ozenith range — 1 to -|- 1 . It 
should be noted here that the energy and angular resolutions between the muons and the neutrinos of the events 
differ significantly for neutrinos and anti-neutrinos due to their different ways of interactions. 

To generate the theoretical data Nj'j for the chi-square analysis, we first generate 500 years un-oscillated data 
for 1 Mton detector From this data we find the energy-angle correlated resolutions (see figs. [T]) in 30 bins of 
energy (in log scale) and 1 bins of cosine of zenith angle ( — 1 to + 1 ) . For a given Ey,we calculate the efficiency 
of having E^> 0.8 GeV (threshold of the detector). For each set of oscillation parameters, we integrate the 
oscillated atmospheric neutrino flux folding the total CC cross section, the exposure time, the target mass, the 
efficiency and the resolution function to obtain the predicted data in the reconstructed logE — cosOzenith grid0. 
We use the CC cross section of Nuance-v3 |i2T| and the Honda flux of 3-dimensional scheme |[28l . 

The atmospheric neutrino flux is not known precisely. There are huge uncertainties in its estimation. We may 
divide them into two categories: I) overall uncertainties (which are flat with respect to energy and zenith angle), 
and II) tilt uncertainties (which are function of energy and/or zenith angle). These have been estimated as the 
following II2TI : 

1. The energy dependence uncertainty which arises due to the uncertainty in spectral indices, can be ex- 
pressed as: 

<I>8,(£)=Oo(£)(£) ' ^ME) 
The uncertainty of 5e =5% and Eq = 2 GeV is considered. 

2. Again, the flux uncertainty as a function of zenith angle can be expressed as 

(cos ej f«<i)o(cosej [1 + 5^1 cos ej.|]. (lO) 

The uncertainty of 8j is considered to be 2%. 

3. A flux normalization uncertainty of 20%. 

4. An over all uncertainty of 10% in neutrino cross section. 

5. An overall 5% uncertainty for this analysis. 

'^One can do this in an another way. This is generating the theoretical data directly for each set of oscillation parameters. To ensure 
that the statistical error is negligible, one needs first to generate a huge number of events. For example, one may generate events for 500 
Mton.year exposure of the detector for each set of oscillation parameters. Then to obtain the theoretical data, one needs to normalize 
the data to 1 Mton.year exposure of the detector dividing the events of each energy and zenith angle bin by 500 since the experimental 
data is considered for 1 Mton.year exposure. This would be the more straightforward method. But the marginalization study with 
this method is almost an undoable job in normal CPU. However, an exactly equivalent result is obtained here using the energy-angle 
correlated resolution function. 



l + 5£logio — 



(9) 
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We consider three types of binning: 

• Type I: The events are binned in the grid of logE — logL plane. We use total number of logE bins ue 
= 30 (0.8 — 50 GeV) and the number of logL bins as a function of of the energy. We consider = 
2x 14, 2x 18, 2x22, 2x26, and 2 x 30 for £ = 0.8 - 1.2, 1.2-2.4, 2.4-3.6, 3.6-4.8, and > 4.8 
GeV, respectively. For the down-going events the binning is done by replacing 'logL' by ' — logL'. The 
factor '2' is to consider both up and down going cases. 

• Type II: The events are binned in the grid of logL — cos Ozenith plane with exactly in the same fashion of 
type I. The only difference is that the binning is done in cosOzenith instead of logL. 

• Type III: The events are binned in 100 log(L/L) bins and replacing 'log(L/L)' by ' — log(L/L)' for 
down-going events. 

For the up-going neutrino, L is the distance traveled by the neutrino from the detector to the source at the 
atmosphere. In case of down-going neutrinos the distance traveled from the source to detector is negligible for 
getting an appreciable oscillation. However, these events help to minimize the systematic uncertainties when 
considered in the analysis. The flux for a fixed E is strongly dependent on the zenith angle. So, for the 
down-going neutrinos, we mapped the zenith angle into L considering the mirror L. This is the same L if the 
neutrino comes from exactly opposite direction. It should be noted here that the angular error makes a much 
smaller error to L when the tracks are near vertical. It increases gradually when the tracks are slanted and very 
rapidly when they are near horizontal. 

For each set of oscillation parameters we calculate the in two stages. First we used t,^ such that ^ = 0, 
which can be obtained solving the equations 11291 . Then we calculate the finalx^ with these values. Finally, 
we find the minimum from these x^ with respect to all oscillation parameters B 



3 Result 



In this section, we first discuss the results qualitatively in a very general way for all analysis techniques. 
Then we compare the results for different techniques. In all cases a global scan is carried out over the three 
oscillation parameters Am^2^ ^23 and 613 for both normal and inverted hierarchies with neutrinos and anti- 
neutrinos separately. We have considered the range of Am^2 = 2.0 — 3.0 x lO^^eV^, 623 = 37° — 54°, and 
613 = 0° — 12.5°. We have fixed other parameters Am2j and 612 at their best-fit values and 5cp = 0. The 2- 
dimensional 68%, 90%, 99% confidence level allowed parameter spaces (APSs) are obtained by considering 

^Here we consider all uncertainties as a function of reconstructed neutrino energy and direction. We assumed that the tilt uncer- 
tainties will not be changed too much due to reconstruction. However, on the other hand, if any tilt uncertainty arises in reconstructed 
neutrino events from the reconstruction method or kinematics of scattering, these are then accommodated in x^- We first incorporate all 
uncertainties in logii — cos O^g^ith bins. Then we re-bin the data in the form what we want, e.g.; logE — logL bins. It should be noted 
that we first binned the data into a large number of cos Ozenith bins compared to number of L bins to get proper binning in logL. 
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Figure 1: The sample energy-angle correlated resolution plots for neutrino (left column) and anti-neutrino (right 
column) for the bins of £v = 0.85 - 0.98 GeV with cos ezenith = -0.40 to -0.20 (upper row) and Ey = 6.84 - 7.86 
GeV and cosQzenith = to 0.20 (lower row). The data are obtained from the simulation of 500 MTon.year 
exposure of ICAL considering no oscillation. 

= Xmin + 2-48, 4.83, 9.43. To obtain the APS in 613 - Am^j i^rnli ~ ^23) plane, we marginalize the over 
023 (613) over its whole range. 

The experiment indicates that the value of 613 is very small compared to 623 |[30l . So, the atmospheric neutrino 
oscillation is mainly governed by two flavor oscillation (v^) \x (Vx). This constrains sin^2623 and |Ani32|. 
Now, there appears a degeneracy in 623 whether it is larger or smaller than 45° due to the sin^ 2623 dependence 
of oscillation probability. However, when the matter effect comes into the play, the effective value of 613 
becomes large and a resonance occurs in (v^) ^ v^. (Vg) oscillation. This breaks the above 623 degeneracy. 

The difference in oscillation probability between two 613 values for neutrinos with NH and for anti-neutrinos 
with IH becomes significant when matter effect comes in the picture (see eq. |3]). We have plotted the APS 
in 013 — Am32 plane considering both neutrinos and anti-neutrinos (i.e. with x^^^^j = Xv +x|) for different sets 
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-0.01 




Figure 2: The oscillation probability of v^, v^. We choose Amjj = -2.5 x 10 ^eV^, 823 = 45° and 813 = 0°. 

of input parameters at 68%, 90% and 99% CL in fig |4l [5] and [6l respectively for each type of binning of the 
data. We see that the matter effect significantly constrains 613 over the present limit, which is a very stimulating 
result for atmospheric neutrino oscillation analysis. 

Again, for the APS in Am^2 ~ ^23 plane, 613 is marginalized over the present allowed range. The APSs are 
shown in fig. |7J[8]and|9]at 68%, 90% and 99% CL, respectively for each type of binning of the data. If the value 
of 613 is nonzero, the matter effect plays a role in determination of the octant of 623 as discussed previously and 
also constrains the 623 range (compare its range for zero and non-zero values of 613). We find that for some 
combinations of (613, 623), the octant determination is possible. 

Now we will compare the APSs coming from different analysis method. From the APSs it is clear that the L/E 
analysis gives very poor results compared to the other two methods. It happens due to the mixing of events 
from different E and L resolutions since the resolution widths are strongly energy dependent. It should be noted 
here that we have not used any selection criteria for the events, which might improve the results. 



8 







log E - log L 
log E - cose,g„j,f, 



































































0.0022 0.0023 0.0024 0.0025 0.0026 0.0027 0.0028 
A (eV^) 



Figure 3: The typical distribution of A^^ witin Amgj. We choose the input of Amgj = +2.5 x 10 ^eV^, 623 = 42° 
and 013 = 7.5°. 



Now we will compare the positive and the negative sides of the rest two methods. We find a relatively stronger 
upper bound of l^m^^ in case of binning in the grids of logS — logL plane than the other case. This is very 
important since it comes from the events with high E and low L values. The L resolution is very poor at low L 
and the statistics is low at high E. However, a stronger bound is obtained for this special type of binning. We 
will explain it with the oscillation probability in vacuum, which is a sinusoidal function ofL/E. For a fixed L, 
the distance between two consecutive peaks in E increases rapidly with E. Again, if we compare the distances 
between two consecutive peaks in E for two fixed values of L, it is larger for smaller L value. Therefore, as one 
goes to smaller L values, this distance in E increases rapidly. So, these two consecutive peaks of the oscillation 
in E can be resolved with much better resolution as one goes gradually from larger L values to lower L values. 
This is pictorially illustrated in fig. |2] To get the reflection of this fact in y}, the finer binning at lower L is 
essential. Though the angular resolution is worsened at lower L, but the rapid increase of E resolution between 
two peaks wins the competition here. This is the main advantage of this type of binning. So, we binned the 
data in a two dimensional grids of logL — \ogE plane §. In type II this behavior is not taken into account in the 
binning of the data. However, there is a disadvantage in type I that the bin size at high L values is very large 
compared to type II, which gives weaker lower bound on Amjj. So, the combination of type I and II (type I at 
the lower range of L and type II for the rest) is a better choice than the individual cases. However, this is not 
studied in this paper, but is reflected when we compare two results. This is also demonstrated in terms of A^^ 
for a typical set of parameters in fig. |3] It should be noted here that the contrast between two methods would 
be prominent when the number of bins in L or cos Szenith will be relatively lowered than that used in this paper. 

For a quantitative assessment of the result, we define the precision P of a parameter t as: 

(J max f iT^in 
^ 
J max _j_ J mm 

^This captures the oscillation effect well in analysis without mixing events from different E and L resolutions. 




9 



Wesee, one can achieve the precision of Am^2 ~ 4.8-7.5% (5.4-8.0%), 6.9-10.9% (8.0-12.6%) and9.6- 
15.7% (10.9 - 17.6%), at 68%, 90% and 99% CL, respectively in type I (II) method. For the input with bi- 
maximal mixing of 623, we find its precision in terms of sin^923 ~ 14.3 — 31.8% (16.9 — 36.9%), 21.6 — 
36.8% (22.4-41.7%), and 28.5 -42.1% (27.9-45.9%) at 68%, 90% and 99% CL, respectively in type I (II) 
method. The precision of 613 is strongly dependent on its input value. For 613 = 0, we find its upper bound 
6.4° (8.0°), 8.0° (9.5°) and 10.1° (11.5°) at 68%, 90% and 99% CL, respectively in type I (II) methods. The 
both lower and upper bounds are also possible for some combinations of (623,613) and it happens mainly for 
623^45°. 



4 Conclusion 



In this paper we have binned the atmospheric data in three ways: i) in the grids of logE — logL plane, ii) in the 
grids of logE — cos6zenith plane, and iii) in the bins of log(L/£'). We have performed a marginalized study 
over Am32, 613 and 623 for neutrinos and anti-neutrinos separately for each method. Finally, we find that in 
spite of very poor resolutions at low L, which is the main problem as Am32 goes to the upper range, one can 
obtain a relatively stronger upper bound in case of binning in logE — logL plane compared to the binning in 
logE — cos6zenith plane. However, it is also found from both analysis that considerable precisions of 613 and 
Am^2 can be achieved and the octant discrimination can also be possible for some combinations of (623,613). 
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Figure 4: The 68% CL allowed regions in Gb Amjj plane for type I (top) type II (middle) and type III (bottom) 
binning of the data with the input of 623 = 40°, 42°, 45°, 48°, 50° with 813 = 0° (first column), 5° (second column), 
7.5° (third column) from neutrinos with NH and 7.5° (fourth column) from anti-neutrinos with IH. 
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Figure 7: The 68% CL allowed regions in 623 - Amjj plane for type I (top) type II (middle) and type III (bottom) 
binning of the data with the input of 823 = 40°, 42°, 45°, 48°, 50° with 813 = 0° (first column), 5° (second column), 
7.5° (third column) from neutrinos with NH and 7.5° (fourth column) from anti-neutrinos with IH. 
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Figure 8: The same plots of fig. Hbut with 90% CL. 
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Figure 9: The same plots of fig. Hbut with 99% CL. 
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